Reduced-symmetry two-dimensional solitons in photonic lattices 
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We demonstrate theoretically and experimentally a novel type of localized beams supported by the 
combined effects of total internal and Bragg reflection in nonlinear two-dimensional square periodic 
structures. Such localized states exhibit strong anisotropy in their mobility properties, being highly 
mobile in one direction and trapped in the other, making them promising candidates for optical 
routing in nonlinear lattices. 
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Nonlinear periodic structures attract increasing atten- 
tion due to the possibility to engineer the transmission 
and reflection properties of waves, opening new applica- 
tions of photonic crystals for all-optical signal processing 
and switching 0] . An important part of such applications 
is the ability to control wave propagation in the form of 
nonlinear localized modes or solitons |2j. Conventional 
approaches are based on wave transport through struc- 
tural defects or embedded waveguides, however increased 
flexibility can be achieved when the internal structure 
and symmetry of the nonlinear state itself select the di- 
rection of propagation in defect-free periodic structures. 

The symmetry of solitons is intrinsically defined by 
the physical mechanisms responsible for light localiza- 
tion, and two distinct scenarios are possible in periodic 
photonic structures with self- focusing nonlinearity. First, 
it is light waveguiding in the regions of high refrac- 
tive index due to the effect of total internal reflection; 
this effect leads to the formation of discrete lattice soli- 
tons 0,0 0,0 0,0- Second, it is resonant Bragg or Laue 
reflection in periodic structures; this effect leads to the 
formation of gap solitons with the frequencies inside the 
photonic bandgaps of one-dimensional |lL IM M. IwL fill. Il2| 
and two-dimensional 0EJE1EI periodic structures. 

In this Letter, we demonstrate that in two-dimensional 
periodic nonlinear systems it is possible to utilize both 
localization mechanisms jointly to obtain self-trapped 
states with different properties along the two principal 
directions in a square lattice. We describe theoretically 
the families of such reduced- symmetry gap solitons and 
study experimentally the nonlinear beam self-trapping. 
Further, we identify the unique mobility of these gap soli- 
tons which is highly anisotropic along the two principal 
axes of the square lattice, suggesting a novel mechanism 
for directional nonlinear wave transport in symmetric lat- 
tices. 

Optically-induced photonic lattices provide an impor- 
tant experimental tool to study the effects of nonlin- 
earity on light propagation in periodic structures, being 



dynamically reconfigurable and possessing strong non- 
linear response at the mW level. We study the prop- 
agation of an optical beam through a two-dimensional 
(2D) optically-induced lattice which is governed by the 
parabolic equation for the slowly varying amplitude of 
the electric field 0,0,111 
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+ f(x,y,I)E = 0, (1) 



where D is the diffraction coefficient and the function T 
describes the refractive-index modulation induced by the 
lattice and the beam. We consider a periodic potential 
imprinted in a biased photorefractive crystal by an inten- 
sity pattern of the form T{x,y, |-E| 2 ) = —"f{h + Ip( x ) + 
|i?| 2 ) -1 . Here It, = 1 is the normalized constant dark ir- 
radiance, I p (x) = I 5 [cos(7r(a; + y)/d) + cos(7r(x — y)/d)] 2 
is the interference pattern [see Fig. ffl L>)] which induces 
a lattice with the period d, 7 is the nonlinear coefficient 
proportional to the applied DC field, and / = \E\ 2 is the 
intensity of the probe beam. 

The propagation of linear (small-amplitude) waves is 
described by the spatially extended eigenmodes in the 
form of Bloch waves ^t|- The Bloch waves are found 
as solutions of the linearized Eq. (JIJ in the form E = 
"4>(x, y) exp(i/3z + iK x x + iK y y), where ip(x,y) has the 
periodicity of the underlying lattice. For a square lat- 
tice the dispersion dependencies (3(K x ,K y ) are periodic 
with respect to translations K x ^ y — > K x ,y ± 2ir/d, and 
therefore are fully defined by their values in the first 
Brillouin zone. Then the calculated bandgap spectrum 
can be folded along the principal axis of the lattice, and 
it is shown in Fig. QJa). The parameters used for the 
calculations are chosen to match the experimental con- 
ditions: the diffraction coefficient is D = z s \/ (4irn x 2 ); 
x s = y s = 1/xm and z s = 1 mm are the scales of the di- 
mensionless variables x, y, z; no — 2.35 is the refractive 
index of a bulk photorefractive crystal, A = 532 nm is 
the laser wavelength in vacuum, 7 = 2.36, and I g = 0.49. 

In a self-focusing medium nonlinearity increases the 
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FIG. 1: (color online) (a,b) Bandgap spectrum and intensity 
distribution of the 2D square lattice. (c,d) Calculated inten- 
sity and phase distribution of the Bloch wave corresponding 
to the X-point of the second band of the spectrum (a). 



beam propagation constant, shifting it inside the gap for 
the modes corresponding to the top of dispersion bands, 
and allowing for the formation of self-trapped states or 
spatial solitons Q- It follows that for a square lattice 
the solitons can originate from the T point of the first 
band and X (or Y) point of the second band 0] [see 
Fig. nja)]. At the T point, the Bloch-wave diffraction 
is fully symmetric along the principal lattice directions 
(x and y), giving rise to discrete lattice solitons [EHI- 
Symmetric superposition of X and Y states gives rise to 
symmetric gap solitons 0, Q| or gap vortices . It 
was predicted 0], that solitons may be also associated 
with pure X (or Y) states, having different widths along 
the principal axis directions. The analysis in Ref. [13( was 
performed using the semi-analytical approach based on 
the coupled-mode theory, which can only provide qual- 
itative description of wide solitons which extend over 
many lattice sites. However, we find that wide solitons 
may experience quasi-collapse and transform into highly- 
localized states. We demonstrate below that their prop- 
erties strongly depend on the position with respect to the 
lattice, and show that this has a fundamental effect on 
the soliton transport. 

The families of gap solitons originating from the X- 
symmetry point are found numerically by the 2D relax- 
ation technique, and they are shown in Fig.^top). There 
exist four different symmetries of the 2D gap solitons, as 
shown in the examples (a-d) for (3 — —2.04. They are 
centered at or between the lattice sites for the horizon- 
tal and vertical axis, respectively. In all the cases, the 
soliton phase is constant in y- and oscillating in the x- 
direction, resembling the structure of the corresponding 
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FIG. 2: (color online) Families of gap solitons shown as nor- 
malized soliton power vs. the propagation constant, (a-d) 
Examples of four types of gap solitons for /3 = —2.04. Lines 
show the position of the lattice maxima. 



Bloch state [see Fig. ^c,d)]. Alternating phase is a key 
signature for localization due to Bragg reflection, whereas 
the constant phase in the y-direction indicates that total 
internal reflection is responsible for self-trapping along 
this axis. We also note that the solitons are broader 
in the x- and more localized in the y-direction, having 
reduced- symmetry with respect to the underlying square 
lattice. This happens because the diffraction is stronger 
in the x-direction due to Bragg scattering, similar to the 
effect of dispersion enhancement in Bragg gratings Q, 
whereas along the y-direction, diffraction is weaker re- 
sembling discrete field tunneling in waveguide arrays . 

The ability to control soliton motion in nonlinear lat- 
tices may open new possibilities for optical switching. It 
was shown that solitons can become trapped in a ID pho- 
tonic structure due to a self-induced Peierls-Nabarro po- 
tential when the optical power exceeds a certain thresh- 
old l2(ij . Most remarkably, we find that reduced- 
symmetry solitons can simultaneously be trapped in one 
direction, and fully mobile in another direction, offer- 
ing a novel approach for control of nonlinear transport 
in symmetric lattices. Indeed, since the powers for the 
types (a) and (b) solitons depicted in Fig. Eftop) coin- 
cide for any value of the propagation constant, it follows 
that the Peierls-Nabarro potential is negligible along the 
x axis, allowing for soliton propagation in this direction. 
On contrary, there is a large difference in the powers for 
the solutions of types (a,b) and (c,d). This implies the 
existence of a significant Peierls-Nabarro potential along 
the y axis and reduced mobility or trapping of the soli- 
tons along this direction. We also note that the states 
(a) and (b) appear energetically preferable and likely to 
be experimentally observable. 

In experiment, we induce a 2D periodic modulation 
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FIG. 3: (color online) Focusing of a beam corresponding to 
the X-point of the second band: (top) width at the crystal 
output in experiment (dots) and theory (lines). Solid line 
and dots - widths along x; Dashed line and triangles - widths 
along y. (a-c) experimental, (d-f) numerical results: (a,d) 
input beam profile; (b,e) output beam profiles at low power; 
(c,f) localized states at high power. Inset - optical lattice. 



of the refractive index in a biased 20 mm long SBN:60 
photorefractive crystal ,16:]. As a light source we use a 
cw frequency-doubled Nd:YVC>4 laser at 532 nm. The 
laser beam is split into two parts by a polarizing beam 
splitter forming ordinary and extraordinary waves with 
respect to the crystalline axes. To create the photonic 
lattice, the ordinary polarized beam passes through a 
diffractive optical element (DOE), splitting it into four 
coherent broad beams. The DOE is imaged by an opti- 
cal telescope on the front face of the crystal, thus creating 
a 2D square interference pattern with a period of 23 fim. 
The lattice extends through the whole length of the crys- 
tal and is oriented at 45 degrees with respect to the bias 
(horizontal) electric field [Fig.|3ftop, inset)], in order to 
reduce the effect of the intrinsic anisotropy of the pho- 
torefractive nonlinearity |2l| . The second extraordinary- 
polarized beam is directed onto a programmable phase 
modulator (1024x768 pixels), to allow for phase and am- 
plitude engineering of the A-point Bloch states from the 
second band. The field structure produced by the modu- 
lator is also imaged onto the front face of the crystal and 
combined with the lattice onto a beam splitter. The front 
and the back faces of the crystal can be imaged onto a 
CCD-camera to capture beam intensity distribution. 



In order to match the field profile of the linear Bloch 
state [Figs.[I{c,d)], we use a broad input beam with sev- 
eral out-of-phase humps [Fig.^a)], positioned in between 
the lattice sites along the x-direction, while centered onto 
a lattice site in the y-direction. To compensate for the 
anisotropic diffraction, the input beam is made elongated 
along the x-direction. Such a broad beam allows for spec- 
trally pure excitation of the Bloch state from the top 
of the second band at the A point, which can then be 
moved adiabatically into the gap as the beam self-focuses, 
thus avoiding radiation into other spectral bands. In 2D 
lattices, the broad beams are often susceptible to col- 
lapse instability, as their width gradually decreases due 
to self-focusing when the input power slightly exceeds the 
threshold for the soliton formation . We find that for 
reduced-symmetry states, the collapse is arrested when 
the beam becomes confined to a single lattice site in the 
y direction, effectively reducing the dimensionality to a 
quasi-lD system where collapse does not occur, similar 
to pulse dynamics in fiber arrays [22I ] . The process of lo- 
calization of the optical beam with increasing of the laser 
power is summarized in Fig. EI The predicted anisotropy 
in diffraction can be clearly seen in Fig. OUtop), where 
the beam widths along the two lattice axes are plotted 
against the input power. The widths are estimated as the 
full width of the half maximum (FWHM) of the Gaus- 
sian envelope of the modulated beam profile in both di- 
rections. At low intensities [20 nW, Fig.[3fb)], the beam 
diffracts linearly and its width increases more in x- than 
in y-direction despite the fact that the input size along 
the x-direction is larger. The spatial structure of the out- 
put beam well represents the characteristic fine-structure 
of the Bloch wave from the A-point of the second band, 
with out-of-phase lobes along the x direction as confirmed 
with intcrfcrometric measurements. At higher laser pow- 
ers, the beam focuses until a strongly localized state is 
formed at about 230 nW [Fig.[3Jc)] which covers approxi- 
mately three lattice sites in the x and one in y directions. 

To support our experimental results and confirm the 
formation of localized gap states we performed direct 
numerical simulations of Eq. with the initial con- 
ditions similar to the experimental ones [cf. Figs. E2a) 
and (d)], defining the input beam envelope as ip(x,y) = 
Aexp(— x 2 /S 2 — y 2 /5 2 )sm(2irx/d), where S x and S y are 
the beam widths along the two lattice axes. We find that 
the results of the corresponding numerical simulations 
[lines in Fig. OJtop) and Figs. EJe,f )] show the same ten- 
dencies as our experiments. Similar to the experiment, 
the width of the output state decreases to a minimum 
value of about 40 /im in the x-direction and down to 
one lattice period (23 /im) in the y-direction. The slight 
increase of the beam width along the x-direction at pow- 
ers higher than 250 nW is due to beam splitting into 
two solitons. This effect is typical for one-dimensional 
solitons [2|, and shows that strongly localized reduced- 
symmetry solitons indeed demonstrate a quasi- ID dy- 
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FIG. 4: (color online) (a-d) Numerically calculated output 
beam profiles for different initial tilts of the beam, (a) no 
tilt, (b) 20 mrad tilt along x, (c) 20 mrad tilt along y, (d) 
20 mrad tilt along x and y. White bars mark lattice periods. 
(e,f) show two experimental examples for the mobility along 
x, the cross marking the beam center at the crystal front face. 

namics. In the numerical simulations, we can also resolve 
the complete beam evolution inside the crystal. The cal- 
culations show that the beam experiences quasi-collapse 
at its initial stage of propagation, whereas it reaches a 
steady-state propagation afterwards. We find that the 
solitons may disintegrate due to the development of weak 
oscillatory instabilities at large propagation distances. 
However, all solutions from the family presented in Fig. El 
propagate steadily for experimentally feasible distances. 

We now investigate the soliton mobility, which is ex- 
pected to be highly anisotropic. We impose an initial 
tilt of the input beam along the different directions and 
study the beam displacement at the output, comparing to 
the immobile soliton position [Fig. |3Ja)]. In simulations, 
an initial tilt of 20 mrad (15% of the Bragg angle) in x 
moves the output for two lattice sites [FigEJb)], whereas 
the same tilt in y leads only to a small deformation but 
no movement of the output state [Fig^Jc)]. Even with 
both tilts superimposed, the soliton moves the same two 
lattice sites along the x axis only [Fig|2[d)]. Two exam- 
ples of the corresponding experimental results presented 
in Fig. Hfe,f) demonstrate the high mobility of the lo- 
calized states along x direction. We underline that the 
lattice itself is uniform in x and y, and the directional 
mobility is only determined by the localized state itself. 
Their ability to robustly move along one row (or column 
for Y state) of the lattice makes them good candidates 
for flexible soliton networks in 2D periodic structures. 
This offers new opportunities compared to the soliton 
networks sugg ested earlier for optical signal routing and 
switching |23| . 

In conclusion, we have presented the experimental ob- 
servation of novel self-trapped gap states in a square two- 
dimensional lattice, where localization is based on the 



combined effects of total internal and Bragg reflection. 
These states originate from the X-symmetry point of the 
lattice bandgap spectrum, and they possess a reduced 
symmetry and highly anisotropic diffraction properties. 
Due to this anisotropy, they exhibit enhanced mobility 
along the modulated direction and are trapped by the 
lattice in the transverse direction, allowing for control of 
directional wave transport with possible applications for 
optical routing and switching in planar nonlinear periodic 
photonic structures. 
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